Reentrant topological phases in Mn-doped HgTe quantum wells 
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Quantum wells of HgTe doped with Mn display the quantum anomalous Hall effect due to the magnetic 
moments of the Mn ions. In the presence of a magnetic field, these magnetic moments induce an effective non- 
linear Zeeman effect, causing a nonmonotonic bending of the Landau levels. As a consequence, the quantized 
(spin) Hall conductivity exhibits a reentrant behavior as one increases the magnetic field. Here, we will discuss 
the appearance of different types of reentrant behavior as a function of Mn concentration, well thickness, and 
temperature, based on the qualitative form of the Landau-level spectrum in an effective four-band model. 

PACS numbers: 73.43.-f, 73.63.Hs, 71.70.-d, 85.75.-d 



I. INTRODUCTION 



The study of topological states of matter has undergone a 
vertiginous growth since the theoretical prediction' "^ and the 
experimental observation"' of the quantum spin Hall (QSH) ef- 
fect. Unlike the quantum Hall (QH) effect, which is generated 
by an external magnetic field, or the quantum anomalous Hall 
(QAH) effect, which requires time-reversal symmetry (TRS) 
to be spontaneously broken without applying an external mag- 
netic field, the QSH state is characterized by TRS and is gen- 
erally driven by the intrinsic spin-orbit (ISO) coupling.^ Nev- 
ertheless, it has been recently shown that in absence of spin- 
flip terms the QSH effect survives even if the TRS is broken. 
This state has been dubbed a weak QSH state^ (or TRS bro- 
ken QSH stated, where the weakness refers to the absence of 
protection by TRS. Indeed, the gap and the topological Chern 
and spin Chern numbers associated with the topological phase 
remain robust if the TRS is broke n by an exchange term^ or by 
an additional magnetic field,'^'^ and there is a quantum phase 
transition to a topologically distinct or to a trivial phase only 
when the gap is closed.!^ 

An interesting open question is what kind of competition 
could originate in systems where there is an externally ap- 
plied magnetic field in addition to intrinsic magnetic mo- 
ments, which on their own would lead to the QAH effect. Re- 
cently, the QAH effect has been studied thoroughly for several 
theoretical models of two-dimensional topological insulators, 
including HgTe quantum wells'" and thin films of Bi2Se;^2]Ull 
doped with transition metal elements such as Mn, Fe, or Cr. 
In addition, graphene has been proposed as a candidate for 
the observation of this effect.^ '"* The influence of a magnetic 
field on Mn-doped HgTe quantum wells has been partially in- 
vestigated in Ref.ilO, with the aim of polarizing the Mn mag- 
netic moment to be eventually able to generate the QAH ef- 
fect upon shutting down the magnetic field. However, here 
we concentrate on aspects not considered so far. Since for a 
certain range of parameters the Zeeman coupling can have a 
similar effect on the Landau-level (LL) spectrum as the ISO 
interaction^ ^ '^ '*> and can also lead to the (TRS broken) QSH 
effect, it is interesting to explore the interplay between the 
usual Zeeman term, which is linear in the applied magnetic 



field, and the non-linear effect arising from the exchange cou- 
pling with the magnetic moment of the Mn atoms. In this 
paper, we show that within a model with a spin-conserving 
Hamiltonian, the TRS broken QSH phase occurs, and that a 
reentrant behavior is present for a certain range of parame- 
ters. 

A reentrant integer QH effect has been experimentally ob- 
served a few years ago in GaAs quantum wells for filling fac- 
tors between v = i and z^ = 4, in the first LL.El The phe- 
nomenon was later understood to occur due to a sequence 
of first-order quantum phase transitions between electron- 
solid (Wigner crystal or bubble phases) and electron-liquid 
phases,1^ and it was grounded on the strong electron-electron 
interactions that dominate the physics at non-integer filling 
factors. Due to the self-similarity of the Hall conductance 
curve,Q2l which displays a fractal behavior, a similar phe- 
nomenon was predicted to occur also for a second-generation 
of composite fermions.^'^ In that case, a series of reentrant 
plateaus would turn out to be quantized at the nearby frac- 
tional Laughlin values.^' 

A second possibility is to observe reentrant integer QH ef- 
fects solely due to the LL structure of the system. For in- 
stance, in Si/SiGe heterostructures, reentrant behavior can be 
driven at the single-particle level by varying the in-plane mag- 
netic field while keeping the perpendicular field fixed, as to 
modify the ratio between the cyclotron energy and the Zee- 
man splitting. ^^ In these systems, the crossings of the LLs for 
spin up and spin down are responsible for the reentrant be- 
havior of the quantized Hall conductivity. For HgTe^^ and 
InAs/GaSb^"* quantum wells, the reentrance of the Hall con- 
ductivity has been used as a practical method to prove the ex- 
istence of such a LL crossing and consequently the inverted 
order of the bands. 

In this paper, we show that by applying a magnetic field 
perpendicular to a Mn-doped HgTe quantum well, the charge 
and spin Hall conductivities may reenter concomitantly, i.e., 
there can be a reentrance of the same topological phase, char- 
acterized by both its charge and spin topological invariants. 
This effect is caused by the nonmonotonic behavior of the 
LL energies due to the nonlinear dependence of the Zee- 
man term on the externally applied magnetic field. A rich 
panoply of LL crossings, combined with the nonmonotonic- 



ity of the LL energies, provides us with regimes of param- 
eters where this reentrant behavior could be experimentally 
accessed. Hgi_yMnyTe quantum wells are ideal candidates 
for the observation of these effects, because they have a strong 
ISO coupling and a large Zeeman g-factor. We use the effec- 
tive four-band model of Re f. [3] to compute the LL spectrurrp2i 
together with the relevant Chern numbers in order to identify 
the QSH state and other QH-like states. Band structure calcu- 
lations are performed for different values of the quantum well 
thickness and of the Mn doping fraction to set realistic param- 
eters for this model. We then study the LL spectra including 
the charge and spin Hall conductivities and determine the con- 
ditions for the observation of reentrant topological phases. 

The outline of this paper is as follows. In Sect. |ll] we de- 
fine the effective model that we use to derive our results. In 
Sect. 



ni 



we compute the LL spectrum, explain the mecha- 
nisms that lead to the reentrant effects, and explore the param- 
eter regimes for which they can be observed. We conclude by 



discussing in Sect. IV the possibilities to resolve the reentrant 



effects in experiments. 



II. THE MODEL 

HgTe and related materials have a zincblende lattice struc- 
ture, so that the physics of the low-energy electro nic ba nds 
is well described by the eight-band Kane model.l^SESI gy 
using perturbation theory (see the Appendix for details), 
higher-energy bands are projected out in order to reduce 
this model to an effective four-band model.'' ''' In this re- 
duced model, the bands under consideration are referred to as 
\Ei+) , \Hi+) , \Ei-) , \Hi-), in this order Here, E and H 
refer to electron- and hole-like bands, respectively, and + and 
— distinguish the two members of each of the two Kramers 
pairs |-Ei±) and \Hi±), hereafter referred to as spin compo- 
nents. The symmetry properties under the parity and time- 
reversal transformations dictate the qu adratic -order Hamilto- 
nian H ^ Hq + Hz + i^cx, witP^ 



fh{k) 
i^ h*{-k) 



(1) 



where 



h(k) = e(k)l2 + djk)a", e(k) = C - D{kl + k^), 

d„(k) = {Ak^,-Aky, M(k)), M(k) = M - B{kl + k^). 

(2) 
Here, (j°' denotes the Pauli matrices and M, A, B, C, and D 
are parameters that depend on the material composition and 
on the thickness of the quantum well. In particular, the vari- 
ations of these parameters induce the topological phase tran- 
sition from a regime where the electronic bands are ordered 
normally to a regime where the order is inverted and where 
the QSH effect is present.^ 

The system is subjected to a perpendicular magnetic field 
Bbz, which we will express in terms of the dimensionless vari- 
able (p, which denotes the magnetic flux per unit cell measured 
in units of the flux quantum h/e. With these definitions, (f> re- 
lates to the magnetic field B and to the magnetic length Zg as 



2n(t> = eBa? /h — a?l^^. For HgTe, with lattice constant 
a — 0.646 nm, the flux value cf) = 10^"^ corresponds to a 
magnetic field of B = 9.91 T. In the remainder of this text, 
we set C = for convenience, and we set a = 1 as the unit 
of length, so that M, A, B, and D all have the dimension of 
energy. 

The materials under consideration show a large Zeeman ef- 
fect, with Lande g factors of the order of 20.^^ We therefore 
consider the Zeeman term in the Hamiltonian, with different 
g factors for electrons and holes. 



Hz = diag(5E,gH, -5e, -5H)(27r0) 



(3) 



where 5e(h) = gE(H)MBVea2 w gE(H) x 9L30meV is a 
rescaled Zeeman parameter, proportional to the Bohr magne- 
ton /iB and to the g factor 5e(h) for electrons (holes). ^-^ 

In quantum wells of HgTe doped with Mn (with molar 
fraction Y, i.e., we consider Hgi_yMnyTe), the presence of 
Mn has a significant effect on the magnetic properties of the 
material. It has been found that for low Mn concentrations 
(Y < 0.07), the material behaves paramagnetically, so that 
its response to the magnetic field is nonlinear: In addition to 
the Zeeman effect (linear in the magnetic field strength), there 
is also a nonlinear contribution from the exchange interaction 
between Mn ions and band states. The exchange interaction 
term is given by?^^ 



iJcx = diag(xE, XH, -XE, -XH)S5/2(Acx27r 



(4) 



where xe and xh are the exchange energies for the electron 
and hole bands, respectively, 

S5/2(x) = I coth(|.T) - i coth(ix) (5) 

is the Brillouin function,'22l 



Aox — 



5 gMnMB^/ea^ 

2 fcB(r + ro) 



5297 K 

T + To 



(6) 



is an exchange parameter, with gMn — 2, and T + Tq is an 
effective temperature, where Tq w 2.6 K.-** Since the elec- 
tron wave function |-Ei±) is a linear combination of the wave 
functions of the F^ and F^ bands, the exchange energy xe is 
a linear combination of the exchange energies Ag and Ap as- 
sociated with these bands, respectively. The only contribution 
to the hole wave function \Hi ±.) co mes from the F^ bands, so 
that xh is proportional to Ap.GSEZJpor a more detailed expla- 
nation, we refer the reader to the Appendix. 

The energy splitting due to the exchange interactions can be 
considered as an effective Zeeman splitting, by virtue of the 
similarity between Hamiltonians Q and Q. Here, one writes 
the Zeeman energy as g^,-^-^2TT<j), where 



5e(h)('/') =.9e(h) 



Xe(h) 

2TT(h 



-B5/2(Aox27r(/)) 



(7) 



is the effective, field-dependent g factor for the electron (hole) 
band. In the low-field limit (27r(/<Acx ^ 1), the effective g 
factor is approximately constant, g'0-^\ (0 — )• 0) = 5e(h) + 
(7/15)xe(h)Acx, derived by using a linear approximation to 



B^/2{x). In the high-field limit 27r(/)Aox ^ 1, the exchange 
interaction energy is almost constant (w Xe(H)) as a function 
of the field, because i?5/2(a;) — > 1 for a: — ?> oo, and as a 
consequence, it depends also very weakly on the temperature. 



III. RESULTS 

In order to derive the LL spectrum, we model the effect 
of the magnetic field Be^ by the Peierls substitution: In the 
Hamiltonian, the momentum /ik is replaced by /ik — eA, 
where A is the gauge potential, such that Bgz = V x A. The 
freedom of the gauge choice allows us to choose the symmet- 
ric gauge, A = (S/2)(— j/,a;,0). Subsequently, we replace 



/C+ — n^.7 



+ iky and fc_ = k^ 



iky by the ladder opera- 



tors a) and a, respectively. These operators raise and lower 
the LL index by 1, and their pre factors are chosen such that 
[a, a)] = iP In the model presented here, we neglect the 
coupling between the two spin bands which would arise in the 
presence of bulk-inversion asymmetry and Rashba spin-orbit 
coupling. By virtue of this decoupling, the two spin bands can 
be treated separately. Thus, the eigenvalues and eigenvectors 
of the Hamiltonian are given by the solutions to the equation 
K{a,a''){\n + l),c\n)) = E^J.l{<l)){\n + I) ,c\n)), with 

the appropriate values for c. Here, the eigenvalues EaL{4>) 
give the energies of the LLs, where n = 0,1,2,... is the LL 
index, cr = +, — refers to the spin components, and i ~ 1,2 
distinguishes between the two solutions that exist for each 
spin component. ^^ ^" For the Hamiltonian that includes the 
(effective) Zeeman effect, the resulting LL energies are given 
hy2MM 



Ej^j^'> = [-2nD -aB+ \a~gf]{2T:(j)) ± ^J[M + {-aD - 2nB + icr.giff)(27r0)]2 + 2nA^{2'K(l)) [n > 1), 
S+,o-M-p + B-.gf)(2^0), E^,o = -M-{D-B + g'^)i2Tr(b), (8) 



where g'^ = g^ ± g^ are the sum and the difference of the 
effective (field-dependent) g factors given by Eq. (|7]i, includ- 
ing the effect of the Mn doping. 

The LL spectra presented here are computed using Eq. ([8]l, 
where the relevant parameters have been derived numerically 
from b and s tructure calculations based on the eight-band Kane 
model,'^5l2ias explained in more detail in the Appendix. These 
parameters have been computed for several values of the 
quantum well thickness d and Mn fraction Y. In particular, 
the dependence of M on F has dramatic consequences: In- 
creasing Y leads to an increase of M, such that it drives the 
system from the inverted regime (M < 0) to the topologically 
trivial regime (M > 0)\^ 

The charge (spin) Hall conductivity in a specific bulk gap 
is defined as the sum of the charge (spin) Chern numbers over 
all occupied LLs below it. Here, by virtue of the decoupling 
of the two spin components in the Hamiltonian, the charge 
and spin Chern numbers of each LL are equal to the sum and 
difference of the Chern numbers C±.„ associated with each of 
the two components. These Chern numbers are well-defined 
due to the spin-conserving nature of the Hamiltonian. Thus, 
the charge and spin Hall conductivity expressed in units of 
their respective quanta, e^/ft- and e/An, are computed as 



cth 



= E(^. 



c_ 



sp _ 



(C- 



+ .n 



C^.n), (9) 



where the summation is over the occupied LLs. These values 
are robust, even in the absence of TRS.** In this model, each 
LL contributes a Chern number of 1, so that the analysis is 
simplified to merely counting the LLs. The presented values 
have been verified by analysis of the edge states in a ribbon 
geometry; see e.g. Refs. l6l and 1761 for the details of this alter- 
native approach. 



r 



The absence of coupling between the two spin states has an 
important consequence for the QSH phase. Since the QSH 
state may be viewed as a combination of two independent 
QH effects for spin up and spin d own, it persists even in the 
absence of time-reversal symmetry.^ -^^ Additional symmetry- 
breaking terms, for instance due to bulk-inversion asymmetry 
and Rashba spin-orbit coupling, would cause an opening of 
a small gap between the edge states, which allows for some 
backscattering in the presence of impurities.!^ 

In Fig. [it a), we have displayed the LL spectrum for an 
undoped (Y — 0) quantum well with d = 7.5 mn. This 
system is in the inverted regime, so that the spectrum dis- 
plays a (weak) QSH gap, (with (cthjCt^) = (0,2)), for 
magnetic fields up to S = 7.6 T, where the LLs cross, at 
(/)^j,ogg = M/[2tt{B - g'f/2)]. In addition to the (weak) 
QSH gap, we observe several spin-filtered (e.g., (cth, cr^) = 
(±1,1)), spin-imbalanced (e.g., {aii,a^) = (3,1)) and or- 
dinary (e.g., {crii,a^) = (2,0)) QH gaps, and a trivial gap 
((cthiO"^) = (0,0)). The (weak) QSH gap is the only gap 
which exhibits a helical edge state structure; all other nontriv- 
ial gaps are chiral. Within this formalism, no other inverted 
gaps form besides the one atip — O, because the involved LLs 
do not cross anywhere else than at = Across- In contrast, a 
tight-binding description of a honeycomb lattice in a perpen- 
dicular magnetic field does allow for other gaps with helical 
edge structures at higher flux and Fermi energy values.* '^ 

The LL spectrum of Fig. flTa) shows two mechanisms that 
lead to reentrant behavior of the Hall conductivity and spin 
Hall conductivity. The first mechanism is illustrated for a 
Fermi energy of 11.0 mcV (lower dashed line), which lies 
just below the energy value at which the two lowest Lan- 
dau levels (LLLs), i.e., the LLs with energies -E+.o and 
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FIG. 1. (Color online) (a) LL spectrum for the quantum well with d — 7.5 nm without Mn doping. Red and blue curves indicate the + and — 
spin components, respectively. For clarity, only the LLs with n < 10 are shown. The numbers inside the gaps indicate the charge and spin Hall 
conductivity, cth and o"^. On the horizontal axes, the flux value (j> and the equivalent magnetic field strength B are given. The dashed lines 
indicate Fermi energies where reentrant QH effect is observed. (b,c) Charge Hall (solid curves) and spin Hall (dashed curves) conductivity as 
a function of the flux (in terms of their respective flux quanta) for the indicated Fermi energies. The red (thick) curves indicate the effect of 
broadening due to disorder (Po — 0.3 meV) and temperature (T = 0.4 K), compared to the results without broadening (black, thin curves). 
(d,e,f) Equivalent plots for d = 7.5 nm and Y = 0.004. The inset in (d) is a magnification of the energy region near the LL crossing. 



E^fi, cross. Holding the Fermi energy fixed and increasing 
the magnetic field, we successively traverse the weak QSH 
gap with ((THjCTy'^) = (0,2), the spin-filtered QH gap with 
(o'HjO'u) ~ (^1' 1)' ^^'^ '■h^ trivial gap, where (cthiO"^) = 
(0,0). Thus, the charge Hall conductivity is for low and 
high magnetic fields, and —1 for intermediate values, which 
characterizes a reentrance of a charge-insulating state (see 
Fig.flTc)). At a Fermi energy slightly above the crossing (e.g., 
E — 13.5 mcV, see Fig.[TJb)), a similar sequence is observed, 
but with a different intermediate state (cth = +!)• In both 
cases, the spin Hall conductivity takes the values 2, 1, and 0, 
and does therefore not show reentrant behavior Clearly the 
reentrance of the Hall conductivity is caused by the structure 
of the spectrum around the crossing of the LLLs. To observe 
the reentrance of the charge Hall conductivity, it is essential 
that the derivatives dE/d(p of the two LLLs at the crossing 
differ in sign, which can happen only in the inverted regime. 
Thus, experimental observation of this type of reentrance pro- 
vides a proof that the HgTe quantum well can indeed be de- 
scribed as an inverted Dirac system."* One may verify that if 
the signs of the derivatives would be equal, then the charge 
Hall conductivity does not reenter Instead, we would observe 
reentrant spin Hall conductivity. We note that crossings of the 
latter type are ubiquitous for higher LLs (n > 0), but they 
are difficult to observe due to the vicinity of other LLs. How- 
ever, later we will show that, under some circumstances, the 
crossings of the LLLs may also be of this type. 

In Fig. [TJd), we show the effect of doping (Y = 0.004) 
on the LL spectrum. Two effects are visible. First, the size 
of the (weak) QSH gap has decreased, consistent with the in- 
crease of M. Secondly, the energies E±fi of the two LLLs 
are no longer linear in the magnetic field. In fact, one of 
these LLLs shows a nonmonotonic dependence on the field. 
As can be observed in Fig. flTd), this nonmonotonic LLL at- 



tains its maximum for a flux value less than Across- Thus, if 
the Fermi energy is located between the energy of the crossing 
and that of the maximum (e.g., if _E = 12.0 meV, see the inset 
of Fig. [ltd) and Fig.flTe)), the spin-filtered QH gap reenters, 
and the intermediate state is the (weak) QSH gap. Thus, the 
system goes from a chiral, to a helical, and back to the (same) 
chiral phase again. This simultaneous reentrance of the charge 
and spin Hall conductivity should be contrasted with the reen- 
trant behavior around the LLL crossing, where only one of 
them reenters, but not both. We remark that such a sequence 
is possible only if the intermediate phase is the (weak) QSH 
phase, and consequently only if the LL involved is one of the 
LLLs, since the higher LLs are all monotonic. Therefore, this 
behavior cannot be observed in the undoped system, where 
the LLL energies are linear 

As can be observed in the inset of Fig. [TJd), the maximum 
of this LLL has an energy close to that of the LL crossing. 
The sequence of charge and spin Hall conductivities is there- 
fore affected by both mechanisms. We shall call this phe- 
nomenon compound reentrant behavior. Above the energy 
of the crossing, the aforementioned sequence (spin-filtered 
QH, weak QSH, spin-filtered QH) is followed by the triv- 
ial gap, so that we get an additional reentrance of the zero 
charge Hall conductivity. Just below the crossing energy (e.g., 
E = 11.3 meV, see Fig.[Tlf)), the sequence of gaps is spin- 
filtered QH (1, 1), weak QSH (0, 2), spin-filtered QH (-1, 1), 
and trivial (0,0). In this situation, the two spin-filtered QH 
phases are different gaps, unlike the sequence above the cross- 
ing. These examples show that the rich compound reentrant 
behavior will appear if the crossing and the maximum of the 
LLLs are close to each other. 

In order to be able to observe the reentrant effects in experi- 
ments, we study the qualitative structure of the LL spectrum as 
a function of the well width d, the doping fraction Y, and the 
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FIG. 2. (Color online) Diagram showing the possible types of reentrant behavior, as a function of the Mn fraction Y and temperature T. In 
the top row, we show diagrams for (a) d = 7.5 nm, (b) d — 6.5 nm, and (c) d — 5.5 nm. The roman numbers correspond to different regimes 
in the qualitative structure of the spectrum constituted by the LLLs, displayed in the bottom row. The characteristics of the regimes (i)-(v) 
are explained in the text. In the top row, the different shades indicate the energy range in which the compound reentrant behavior is present, 
compared to the LL broadening: Brighter colors indicate a larger range, making the effect easier to observe. The dashed line indicates the 
parameter values where the crossing and the maximum of the nonmonotonic LLL coincide, separating (ii) and (iii). The shading in the bottom 
row indicates the weak QSH (red/dark gray), spin-filtered QH (blue/light gray), and trivial (white) gaps. 



temperature T. More specifically, for a fixed choice of param- 
eters, we analyze whether one of the LLL is nonmonotonic, 
and whether the LLLs cross. Furthermore, if the nonmono- 
tonicity and crossing appear at the same time, we determine 
the position of the maximum/minimum and the crossing rela- 
tive to each other. For simplicity, we restrict ourselves to the 
structure of the two LLLs. 



The bottom row of Fig.l2]displays five qualitatively differ- 
ent LLL spectra, which distinguish the regimes as given by 
Fig. l2ta)-(c). These regimes are characterized as follows. 
For regimes (i)-(iii), the band gap has inverted order (i.e., 
M < 0), and therefore shows the QSH phase at zero mag- 
netic field. In regime (i), the LLLs are monotonic, so that the 
only mechanism that leads to reentrant effects is the cross- 
ing. In regimes (ii) and (iii) one LLL is nonmonotonic, so that 
we have compound reentrant behavior. These two regimes 
are distinguished by the flux value of the maximum, which is 
smaller (ii) or greater (iii) than the flux value of the crossing. 
In the case (iii), both LLLs are increasing at the crossing, so 
that we observe reentrant spin Hall conductivity, as argued be- 
fore. In regimes (iv) and (v), the band gap is normally ordered 
(i.e., M > 0), so that we find a trivial phase at zero magnetic 
field. In this situation, the LLLs do not cross, and the only 
mechanism that can lead to reentrant behavior is the presence 
of a nonmonotonic LLL, as is the case (iv). For regime (v), 
both LLL are monotonic and do not cross, thus preventing 
any type of reentrant behavior. 



IV. DISCUSSION 

Let us finally comment on the ability to resolve these reen- 
trant effects, based on the range of the Fermi energies for 
which they are present. In order to estimate the observabil- 
ity, we compare this energy range to the broadening of the 
LLs, that will cause the change of conductivity across a LL to 
be smooth rather than step-like. Here, we consider a gaus- 
sian broadening with width F,^* that incorporates both LL 
broadening due to disorder and the smooth variation of the 
fermionic filling function at finite temperatures. The broad- 
ening due to disorder has a field-dependent width Fdis = 
Tq^B/Bo, where Fq - 0.1-2 me V and Bq = 1 T.l^^ The 
thermal broadening is approximated by a gaussian with width 
Fth = \/2/37rk-QT. InFig.llb,c,e,f), we illustrate the smooth 
transitions of the conductivities due to the combined effects of 
both types of broadening. 

We consider the compound reentrant effects displayed in 
Fig. |2|ii) and (iii) to be well resolvable if the difference be- 
tween the energies at the crossing and at the maximum ex- 
ceeds twice the broadening width 2r. Indeed, in that case, 
the difference between the actual conductivity values and the 
quantized ones in absence of broadening effects is < 0.08. 
However, the value of 2r is not a hard limit: Variations of 
the quantized Hall conductivity, even if they are far from the 
quantized values, may already be considered as a signature for 
a reentrant effect, for example as demonstrated in Fig.[T|e). In 
the diagrams of Fig. l2|a)-(c), the different shadings indicate 
this distance compared to F. In the brightest regions, the dis- 



tance between the LLs is larger than 2r, sufficient for the reen- 
trant effect to be observed. We find that for Fq = 0.3 meV, 
situation (ii) is difficuh to observe due to the small energy 
difference between the LLs, whereas the observation of (iii) 
is easier close to the critical doping, above which the system 
goes to the trivial regime (iv). For the observation of the com- 
pound reentrant effects, thicker wells are favorable because 
the energy range where the effects appear is larger. The simple 
reentrant effect due to nonmonotonicity, as in situation (iv), is 
generally present in a large energy regime and therefore its 
observation is less affected by the LL broadening. 

Transport experiments with HgTe quantum wells have so 
far concentrated on the charge Hall conductivity of the system. 
For example, the reentrance of the charge Hall conductivity 
has been utilized to identify the possible regime of the QSH 
phase. ^^ Observation of the simultaneous reentrant behavior 
of the charge and spin Hall conductivity would also require the 
availability of a spin-sensitive detector, e.g., a contact consist- 
ing of a tunneling barrier and a ferromagnet.^^ However, this 
technique has the drawback that it only works at low fields, 
within the hysteresis range of the ferromagnet. Another detec- 
tion mechanism could be a local magneto-optical Kerr effect 
(MOKE) experiment,^^ although this measurement would be 
difficult due to the small bandgap of the semiconductor. The 
inverse spin Hall effect may provide a way to measure the 
spin polarization of the edge states. ■''* Nevertheless, the mea- 
surement of the charge Hall conductivity at multiple Fermi 
energies together with knowledge of the structure of the spec- 
trum may provide indirect evidence for the existence of these 
reentrant effects. 

In conclusion, we have demonstrated that the nonmono- 
tonic behavior of the LLs in the presence of Mn doping leads 
to reentrant topological phases, and that the vicinity of LL 
crossings leads to rich compound reentrant behavior. Five 
different qualitative forms of the structure of the LLLs were 
shown to occur in the parameter space characterized by Mn 
doping, well thickness, and temperature. Furthermore, we 
have investigated the effects of LL broadening to estimate the 
ability to resolve the reentrant effects in experiments. 
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Appendix: Numerical methods and derivation of the four-band 
effective model by perturbation theory 

In this appendix, we illustrate the used numerical method 
and relate it to the perturbation theory which allows us to de- 



termine the parameters of a four-band effective model. In the 
Kane model, the band structure of the material consists of 
eight bands. ^S However, the two bands jF^, ±1/2) are sep- 
arated by approximately 1 eV from the other six bands and 
will be neglected here. The resulting six-band modified Kane 
Hamiltonian is written in the basis jF'^, 1/2), |r'',-l/2), 
|^^3/2), |r8,l/2), |r8,-l/2) and |r8,-3/2), which we 
denote as |1), |2), |3), |4), |5) and |6) for short in the follow- 
ing. The Hamiltonian can then be written as 

H = Ho + H2 + H,^, (A.l) 

where Hq is the six-band Kane Hamiltonian,'2512SI H^ is the 
linear Zeeman term and iJox is due to the exchange interaction 
between the Mn ions and the band states in a magnetic field B 
in the z direction. The Zeeman term reads 



H. 



9o „ 
z == yMBOCTz 



HI = k'^ibBJ,, 



(A.l) 



for the (decoupled) conduction (|1) and |2)) and valence (|3), 
|4), |5), and |6)) band parts of the Hamiltonian. Here Jz is the 
angular momentum operator, go is the bare Zeeman g-factor 
of HgTe, and k' is a phenomenological parameter. ^^ The ex- 
change term, induced by the sp-d coupling between the Mn d 
level electrons and conduction or valence band electrons, has 
a similar form as the Zeeman term, and reads 



m. = -A, 



^cx - ^^\Jz 



(A.3) 



where A, = 0.2 eV x Y{S) and A^ = -0.3 eV x Y{S) 
are the coupling constants between the Mn spin S and the 
conduction band (A^) or the valence band (Ap), respectively, 
and Y is the mole fraction of Mn^+ ions. The polarization of 
the Mn spin S is assumed to be in the z direction. We regard 
the Mn spin as a classical spin and use the mean field value 
(S) instead of S, which yields 



S={S)= -5oB 



5/2 



5gMnfJ-BB 

2kBiT + To] 



(A.4) 



where i?5/2 is the Brillioun function as given by Eq. (|5j, Sq — 
5/2, gun = 2, and Tq w 2.6 K for Mn.lSllThe argument of 
,65/2 in this equation is equal to Aox27r(/), cf. Eqs. Q and (|6]l. 
Now, we consider the above model in a periodic superlattice 
grown in the z-direction with well width d and barrier width 
L — d. In the limit of large L — d, it becomes equivalent to a 
single quantum well. Due to the periodic boundary condition 
along the z-direction, according to Bloch's theorem, we can 
write the wave function as 



27r 



(A.5) 



where k = (k||,/c^) = ik^,,ky,kz) and (r||,z) = {x,y,z). 
The in-plane wave vector ky is a good quantum number for 
the system, and k^ is the superlattice wave number in the z di- 
rection, taken to be zero, because the quantum wells are effec- 
tively decoupled for large barrier thickness L — d. U^^{z) is a 



d = 7.5nm, y = 0.02 



> 



cq 




FIG. 3. (Color online) Comparison of tlie energy dispersion calcu- 
lated from the full Hamiltonian (blue solid curves) and from the ef- 
fective model (red dashed curves) for d — 7.5 nm and Y — 0.02. 
The two results exhibit a good overlap in the low-energy regime, 
which demonstrates the reliability of the effective model. 



multi-component periodic wave function U^{z + L) = U^{z) 
of the ^-band, which is expanded in terms of a plane-wave 

basis as 



|c^^W> = E<aKA) = 



n.A 



E< 

n,X 



1 



^i{27rn/ L)z 



12-K 



|A>, 



where |A) denotes the component A = 1, 



(A.6) 
, , 6 of the wave 



function, and the expansion coefficients a^^ ^ are functions of 

kj|. The eigenequation for these states is given by H^^ = 
E^'i^, where E^ depends on k||. With the expansion ( |A.6| l, 
we find 



E 



{n,X\H\n',y)ai,,, 



Efd 



?"«,A- 



(A.7) 



A truncation method is applied and a finite number of basis 
vectors {n ~ —N, —N + 1, . . . ,N — 1, N) is used to solve 
this eigenvalue problem to obtain the coefficients a^ ^ . Given 
the fact that we are only interested in the low-energy physics, 
taking A^ = 20 yields a solution that is sufficiently accurate. 

Next, we relate the perturbation theory to the previous nu- 
merical method. The Hamiltonian dA.lb is divided into 



H — iJk„=n 



H, 



(1) 



(A.8) 



where iJk|i=o is treated as the zero-order Hamiltonian and 
H^ as the perturbation. The wave function at the T point 
(k|| = 0) can be obtained from the numerical calculation, 
which is denoted as 



|c^LoW> = E/«-^wi^) 



(A.9) 



with f^.x{z) — 'J2n'4i \\^)- Only the subbands 
\Ei+),\Hi+),\Ei-),\Hi-), which are denoted as 
|A), |i?), |C), |D) for short, are concerned in this calcu- 
lation. Using symmetry arguments, we obtain 



\A)-- 


= /^.i(z)|1)+/a,4(z)|4), 


\B)- 


= /b,3W|3), 


\C)-- 


= /c,2W|2)+/c,5W|5), 


\D)- 


(A. 10) 



where /a,i = /c,2,/a.4 = fc,5jB,3 = fos- Under 
two-dimensional spatial reflection, Ja,!, fc,2, /s,3: /d,6 have 
even parity and /a, 4, fc'.5 have odd parity. Furthermore, in or- 
der to take into account the contribution of the other subbands 
in second-order perturbation theory, additional states \E2±), 
\LH±), \HH2±), and \HH3±) (the second electron, light 
hole, and second and third heavy hole bands, respectively) are 
also solved numerically and can be written in a similar way.l^ 
With the obtained zero-order wave function, we apply the 
second-order perturbation formalism^^ 



i/|m)+^i(m'|iJ«|s)(.|i/(;V) 



-Ljf-n/ -LJs J^m -^s 



(A.ll) 



to obtain the effective model given by Eqs. ([TJ-Q. Here 
|m), \m') are the states chosen from \A), \B), \C), and \D) 
while \s) is one of the intermediate states |i?2±), \LH±.), 
\HH2±), and IHH^i). With this approach, we relate the 
parameters of the effective model Q-Qito the parameters 
of the six-band modified Kane model ( |A. l| l-( |A3| l. We find 
that for the effective mass parameters D and B and for the 
effective g factor ^e, we need to take into account the second- 
order perturbation, while for the other parameters the first- 
order term is accurate enough for our purpose. As the deriva- 
tion is straightforward and the expressions for the parameters 
are quite lengthy, we do not write them explicitly here. As an 
example, the exchange parameters xe and xh are given by 



XES5/2(Aex2^0) = -(FiA, 

XHS5/2(Aox27r0) = -A^ 



FiAp/i), 



■'pi 



(A. 12) 



where Fi = {fA^\fA,i) and F4 = {fA,4\fAA), and Acx27r0 
is the argument of ^5/2 in Eqs. (|4]) and ( |A.4| i. In Table II] we 
show the numerical values of these parameters for severaT dif- 
ferent well thicknesses and different Mn doping. In Fig.p] the 
energy dispersion calculated from the effective model using 
the parameters in Table |l] is shown to fit well with that calcu- 
lated from the full Kane model at small k. This result justi- 
fies the use of the effective model to discuss the low-energy 
physics, in particular in the energy range where the reentrant 
behavior occurs. In this paper, we have restricted ourselves 
to wells with a thickness d < 8.1 mii, because above this 
value, the H2 band lies between the Ei and Hi bands, and 
in that case the four-band model is no longer accurate, espe- 
cially in the energy regime of the valence band. Nevertheless, 
the mechanisms for appearance of the reentrant effects may 
still be present for thicker wells. 



TABLE I. Parameters for the four-band effective model [Eqs. (TJ-l|4}], obtained by perturbation theory 


from the full Kane model 




d (nm) 


Y 


C (meV) 


M (meV) 


A{eV) 


B(eV) 


D(eV) 


ffE 


Sh 


Fi 


^4 


5.5 


0.00 


-16.9 


8.8 


0.60 


-1.15 


-0.73 


15.8 


1.22 


0.62 


0.37 


5.5 


0.01 


-5.8 


20.0 


0.62 


-1.05 


-0.63 


14.4 


1.29 


0.64 


0.35 


5.5 


0.02 


5.6 


31.5 


0.64 


-0.96 


-0.55 


13.2 


1.36 


0.66 


0.33 


5.5 


0.03 


17.4 


43.3 


0.66 


-0.89 


-0.48 


12.2 


1.42 


0.68 


0.31 


6.5 


0.00 


-24.4 


-4.9 


0.58 


-1.45 


-1.04 


20.0 


1.22 


0.58 


0.41 


6.5 


0.01 


-13.9 


5.7 


0.60 


-1.30 


-0.88 


18.0 


1.28 


0.61 


0.38 


6.5 


0.02 


-3.0 


16.6 


0.62 


-1.17 


-0.75 


16.3 


1.35 


0.63 


0.36 


6.5 


0.03 


8.4 


28.0 


0.65 


-1.06 


-0.65 


14.9 


1.42 


0.66 


0.34 


7.5 


0.00 


-29.9 


-14.6 


0.55 


-1.87 


-1.45 


24.3 


1.21 


0.55 


0.44 


7.5 


0.01 


-19.9 


-4.6 


0.58 


-1.62 


-1.20 


21.8 


1.28 


0.58 


0.41 


7.5 


0.02 


-9.5 


5.8 


0.61 


-1.42 


-1.00 


19.5 


1.34 


0.61 


0.39 


7.5 


0.03 


1.4 


16.8 


0.63 


-1.26 


-0.85 


17.6 


1.41 


0.64 


0.36 
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